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Abstract 

It has been widely accepted that relativistic effects 
have to be used for modeUing the wide classes of 
events associated with active galactic nuclei and 
galactic super luminal sources. Here we propose a 
simple analytical model of decelerated motion which 
is able to explain the shape of double-peaked light 
curves of 7-ray bursts. We show that the differences 
between individual bursts of this class may be due to 
the angle between the line-of-sight and direction of 
the outflow, and that the duration of the burst is in- 
timately connected with the deceleration parameter. 
Assuming special relativity formulae and the orien- 
tation of the observer with respect to the outflow, 
we are able to model observed double-peaked light 
curves and obtain good quality fits. We propose 
that such smooth double peaked bursts could be pro- 
duced by new born galactic super luminal sources. 
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1 Introduction 

Only a small fraction of 7-bursts detected by 
BATSE (onboard the Compton Gamma Ray Obser- 
vatory) possess an optical identification. Therefore, 
most of the registered bursts do not have predictions 
for their total energy and we can expect that among 
unidentified bursts may exist some which happened 
close to us. Many active galactic nuclei (AGN) ex- 
hibit relativistic jets, and the Galactic jet sources, 
or microquasars, also exhibit outflow in the form 
of jets (Mirabel & Rodriguez IQ ). We can expect 
that lower energetic type of 7-ray bursts, may be 
connected with a class of new born galactic super- 
luminal sources as it may by inferred from a recent 
review by Ghisellini and Celotti (0] see also \V2.\). 
However, regardless of whether the radiation from a 
burst is isotropic or anisotropic, an observer would 
see only a small fraction of the radiation emitted in 
the burst, even if the observer is looking directly at 
the beam which has finite angular width, he still sees 
only a fraction of the beam. The majority of radi- 



ation detected by the observer will be that having 
the maximum relativistic intensity boosting in the 
direction of the observer. Thus, the shape of the 
burst afterglow light curve should depend on beam- 
ing effects (Rhoads _17_). 

The Lorentz time dilation between the frame of 
reference of the burst and that of the observer, mod- 
ifies the time-evolution of the observed flux of pho- 
tons. It is known that radiating plasma eventualy 
decelerates (Piran Such a situation requires 

that we should derive a simple and fully analytical 
model, which allows us to investigate the effects pro- 
duced in light curves by deceleration. We show be- 
low that the simplest model having constant decel- 
eration can explain particular light curves obtained 
with BATSE. The bursts detected by BATSE con- 
sist of either single, double or multiple peaks in the 
light curve, but we should realize that such a classi- 
fication is determined by our resolution capabilities 
in the time domain. Because we are interested in 
the influence of deceleration we will avoid discussing 
any spatial structure of the source and we will as- 
sume that the temporal resolution of observational 
data is does not allow us to detect spatial extension 
of the emitting region. In the following, we show 
that light curves exhibiting double peaks can be ex- 
plained in the framework of a simple model in a sim- 
ple way. The double peaked light curve appears to 
be a generic property of models with a decelerated 
motion. 

2 The model and its properties 

Consider the motion of a single blob of plasma in a 
relativistic jet emerging from an active region (c.f. 
Dar & Riijula "T). Suppose that this has the sim- 
plest form of decelerated motion, i.e. motion with 
constant deceleration 5 in a straight line - the simpli- 
est assumption that could be made. The agreement 
we find with observations (see Sec. 3), supports this 
assumption although the physical reason for con- 
stant deceleration is not clear. One-dimensional an- 
alytical calculations give us (Rybicki & Lightman 
TS ) the time dependence of the Lorentz factor 7 



and velocity (3 — v/c in the observer's frame: 
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where 70 = jit = 0), and /3o = P{t = 0). We also 
use the Doppler correction of light travel times in 
the form: 
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where 6 is the angle between direction of motion 
and the Hne of sight, t — (70/^0 — lP)c/g, and r — 
(70 — 7)c^/5) are quantities in the observer's frame, 
and Tabs is the time measured at the location of the 
observer. 

The second assumption made is that the tempo- 
ral evolution of the flux of photons emitted during 
the burst is described by a power law F ^ j/^^T^^ 
in the rest frame of the blob. Such a power law de- 
scription of flux evolution is widely discussed and 
frequently (assumed) obtained as an approximation 
of various emission mechanism in the literature (van 
Paradijs et al. Therefore the flux of photons 

received by the observer is: 
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where the time elapsed in the frame of the emitter 
T and the Doppler boosting factor D are given by: 
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The shape of the light curve implied by Q is 
determined by the values of cos0, a, /i, and 79. 
The deceleration parameter (g/c) has no influence 
on the curve shape but only on the burst duration. 
For any set of these parameters we can simply cal- 
culate that there is a T^jJ°^ after which the blob 
stopped its motion. From (2), P — imphes that 
^o™"^ = [7o/3o-(7o-l) cos6']c/5, providing the value 
of TX'"'. When the blob has stopped (Fobs > T^T) 
we adopt to describe the burst simply by putting 
D = 1. At this point, times intervals measured in 
the source and observer frames are equal. Thus we 
are able to model light curves in terms of the pa- 
rameters: cos 6*, a, /i and 79. 

Figure 1 shows the properties of our model and 
the influence of various parameters on the shape of 
the light curve, including the angle between the ve- 
locity vector of the blob and the line of sight. In the 
case of small angles, we obtain a single burst with 
a 'broken power law' shape of light curve in the ob- 
server's frame, studied by several authors (e.g. Sari 



Table 1: Fitted values of parameters. In paretheses 
are values obtained by independent fitting of the 
model with the asumption ^end = 10 (see text). 



GRB 


990531 




000519 




BATSE trigger 


7975 




8111 




70 


15.2 


(119.) 


7.0 


(35.8) 


cos6' 


0.183 


(.985) 


0.070 


(.847) 


a 


1.42 


(1.97) 


-0.15 


(1.36) 




3.11 


(3.18) 


1.46 


(1.65) 


g/c 


0.489 


(.067) 


1.113 


(.747) 


To 


2.071 


(2.02) 


0.679 


(.611) 



at al. [201, Rhoads ^21, Panaitescu & Kumar [TK]1. 
For larger angles we obtain no break in the light 
curve slope but simply a brightening due to the de- 
beaming effect (of the radiation) at later times in the 
decay leading to the formation of a second observed 
peak. Thus, a double-peaked burst can be produced 
by a single blob. The value of critical angle demar- 
cating these effects depends on the spectral index a 
of the burst and its decay index /x. One can also see 
that a larger 70 leads a to a larger ratio between the 
second peak intensity and the minimum intensity 
between the two peaks. 



3 Two examples of fitting 

To test the model, we selected two double-peaked 
7-ray light curves from the BATSE archive (Mal- 
lozii & Six j7|) which are thus consistent with be- 
ing produced by a single blob. We chose the bursts 
which took place on 1999, May 31 and 2000, May 
19 designated GRB 990531 and GRB 000519. We 
fitted the data in each case with the above model 
defined by iQ} using a fitting procedure in which 
each of the 4 free parameters is assigned random 
values by a Monte Carlo procedure, this being re- 
peated many times until an acceptable fit is obtained 
(i.e. RMSfit < 3 X RMS of the fluctuation be- 
fore the burst). Because we adopted the simplest 
model of evolution with a singularity at r = 0, (i.e. 
Tabs = 0, - the trigger time) we were unable to cal- 
culate the growth phase of the burst and we fitted 
only its descending phase from the first peak (from 
Tabs = Tq > 0) and the entire second peak. The 
results of fitting are shown in Figure 2 & 3 and in 
Table 1. 

In both events, the shape of the light curve near 
the second peak (see Figs. 2 & 3) is determined 
physically by the debeaming process (as explained 
above in the discussion of Fig.l) when most of the 
bulk kinetic energy was transformed into thermal 
energy of the surrounding matter via shock interac- 
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Figure 1: Properties of the model i.e. the depen- 
dence of a hght curve shape on parameters /i, F, 
cos 9, and a. Heavy curves are those with fi — 1.2, 
a — 0.7, 7o — 7.2, cos 9 = 0.5. The curves end at 
T = r"(J°^ i.e. at the time when the blob stopped 
its motion. For different values of 70 and cos 9 we 
change the value of g/c in order to make the same 
value of TJJJ"^ for each curve. 



tions, according to the model of Meszaros & Rees 
( 9 , see also Piran We suggest that this is 

the time when optical and/or radio emission could 
rise to a detectable level. Because in the case of 
strong deceleration, the second peak is located a 
few seconds after the first peak, we are unable to 
see their afterglow as optical/radio transients since 
there is not enough time after the high energy de- 
tection to move the telescopes to the correct coordi- 
nates (Paczyriski & Rhoads ^^). 

4 Discussion 

Our model describes a point-like element of plasma 
moving with relativistic speed. Using a Lagrangian 
description of fluid motion we could construct a 
model of more complicated flow similar to a super- 
nova shell during its explosion (Colgate 1 , Galama 
et al. ,3 ). In this case we would observe differ- 
ent parts of the shell possessing different values of 
the parameter cos (6') and the observed light curve 
shape would be a superposition of the single burst 
light curve models. Since the light curve data are flt- 
tcd well using a model involving only a single blob 
of plasma, we do not have to assume that the shape 
of the moving volume influences the observed light 
curve as proposed by Moderski et al. ( 11 ). 

We realize that the low value of the Lorcntz fac- 
tors derived in this paper are inconsistent with all 
but the very softest gamma-ray spectra (Woods & 
Loeb [22, Lithwick & Sari (i ). However, our model 
may be interpreted from the viewpoint of an in- 
ternal shock model scenario (Kobayashi at al. 0, 
Piran jj^li Dar & Rujula |2j) which assumes that 
moving blobs (or shells) possess different velocities 
which lead to collisions between them. Each col- 
lision generates a relativistic shock which radiates 
and decelerates continually. In order to compare 
this scenario with our model we need to introduce 
a new parameter 7end (which is the Lorentz factor 
after the deceleration is switched off), but we are 
unable to determine its value because the correc- 
tions produced by 7end take the form of a scaling 
factor (relative to g). This means that we are able 
to determine it only when we have some predictions 
of total (peak) power or the true time scale of the 
burst phenomenon. In the case of a burst having an 
afterglow, it is obvious that 'jend has a value much 
greater than 1. Such a value may be calculated as 
the initial value in the model of a particular after- 
glow. Therefore values given in Table 1 should be 
scaled in the case of comparison with true physical 
conditions inside the shock and 70 represents the 
difference in Lorentz factors between two collided 
blobs. In particular, this should be done in the ma- 
jority of GRBs which could not be explained by a 
low value of 70. As an example, in Table 1 we put in 
parentheses the values obtained during independent 
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Figure 2: The observed and the best fitted light 
curves for 7-ray burst T7592. The model has param- 
eters: 70 = 15.2, ^ = 3.11, a = 1.42, cos(6l) = .183. 

fitting of the model in the case of ^end = 10. 

Given that GRBs will have a distribution in space, 
it is clear that some bursts will take place relatively 
close to us. If we accept that the bursters produce an 
extreme relativistic outflow, it is more likely that the 
outflows will move in a direction different from the 
line of sight to the observer, compared with a more 
distant burst which we only see if there is alignement 
(Paczyhski 13 J. In such a case, we will see the burst 
light curves that are strongly modified by relativistic 
effects, but with small boosting in the luminosity. 
Because we obtain small values of cos 9 and 70 by 
fitting, we deduce that the two examples of light 
curves presented in Figures 2 & 3 belong to this type 
of burst and they are much closer to us than bursts 
for which redshifts have been determined (Meszaros 

Since we obtain low values of cos 9, therefore the 
outflow direction is far from the line of sight and 
we have no Doppler boosting of radiation in the ob- 
server's frame. Energy considerations lead us to the 
conclusion that the examples of light curves discused 
here, were produced by objects close to us than those 
detected at cosmological distance. If we have a burst 
event at a distance of a few kpc in comparison to 
Gpc, the flux can be (10^^ times larger which allows 
us to see an event which is not boosted but dimmed 
by relativistic effects (or may be also less energetic). 
A scaling factor of 10^^ allows us to calculate the 
upper limit for the Lorentz factor as: lO'^^ > {7^)^ ; 
i.e. 7 < 10^'^ ~ 30. We use 7^^ twice because firstly 
we have no Doppler hosting and secondly because 
we have Doppler dimming. 

It is well known that galactic microquasars (like 
GRS1915-I-105) eject blobs with relativistic speed 
(Mirabel & Rodriguez ^3)- If we assume that sim- 
ilar sources are born in a much more dense environ- 



Figure 3: The observed and the best fitted model 
light curves in 7-ray burst T8111. The model has 
parameters: 70 — 7.0, /i = 1.46, a = —0.15, 
cos(e') = .070. 

ment and work in a similar manner ejecting blobs 
which are abruptly stopped leading to an explosive 
transfer of kinetic energy into radiation, we have a 
physical interpretation of the results obtained. 

5 Summary 

This paper discusses the simplest model of the in- 
fluence of deceleration on the observed shape of the 
light curves of 7-ray bursts. Taking into account 
the deceleration, we obtain important modification 
of the power law decaying light curve. The trans- 
formation of times between the proper frame of a 
blob and the observer's frame, i.e. a logarithmic 
dependence (Eq. 5), may be interpreted as an am- 
plification of the time interval, i.e. at the begining 
of the burst we can observe processes lasting few 
milliseconds as events lasting a few seconds. 

In any case of emission, a spatialy extended region 
could be modelled as a superposition of point-like 
sources. Therefore the shape of the light curve of 
a point-like source give us the limit for the fastest 
temporal variations of the flux. This fact will lead 
us to the prediction that all spatially extended flows 
should possess more smooth light curves than pre- 
dicted by our model. In the framework of a deceler- 
ated motion model, we are able to produce shorter 
bursts (i.e. with lower value of TJJJ"^) if we assume 
stronger deceleration. All events shorter than the 
pulse of radiation discussed here, i.e. induced by 
a deceleration process, probably come from single 
events of coherent emission, but for discussion of 
such an event we should first separate the influence 
of special relativity effects from observational data. 
Therefore when one postulates a spatial shape of the 



outflow (e.g. jet type), the temporal properties of 
the emiting region should be discussed first with in- 
clusion that the decelerated motion gives us the pos- 
sibility of seeing the millisecond events expanded to 
more than a few seconds. 

Values of 70 obtained by fitting suggest that the 
events under consideration belong to a lower en- 
ergetic type of 7-ray burst, which may be con- 
nected with a class of new born galactic superlu- 
minal sources (c.f. |HIT^). 
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